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ABSTRACT 



Context. Rotation curves of interacting galaxies often show that velocities are either rising or falling in the direction of 

the companion galaxy. 

Aims. We seek to reproduce and analyse these features in the rotation curves of simulated equal-mass galaxies suffering 

a one-to-one encounter, as possible indicators of close encounters. 

Methods. Using simulations of major mergers in 3D, we study the time evolution of these asymmetries in a pair of 

galaxies, during the first passage. 

Results. Our main results are: (a) the rotation curve asymmetries appear right at pericentre of the first passage, (b) the 

significant disturbed rotation velocities occur within a small time interval, of ~ 0.5Gyr h~^, and therefore the presence 

of bifurcation in the velocity curve could be used as an indicator of the pericentre occurrence. These results are in 

qualitative agreement with previous findings for minor mergers and fly-byes. 
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1. Introduction 

Rotation curves of interacting and merging galax- 
ies are often highly disturbed. First observations of 



I Rubin et al. (1991)| showed several cases of sinusoidal ro- 
tation curves and asymmetries for a large fraction of com- 
pact gro up galaxies obtained from long-slit spectroscopy. 
Although Mendes de Oliveira et al. (2003)| have shown that 
many of these peculiarities are smoothed when rotation 
curves (hereafter RC) are derived from 2D velocity maps, 
some disturbances still remain, indicating their true nature 
as a global feature of the galaxy. For interacting galaxies, 
features on the RCs are related not only to the structure of 
the disk itself, but also to the external perturbations due 
to the companion (e.g. Rubin and Ford 1983; Elmegreen 
and Elmegreen 1990; Chengalur et al. 1994; Salo andLau- 
rikainen 2000; Marquez et al. 2002; Hernandez Toledo et al. 
2003; Fuentes Carrera et al. 2004; Garrido et al. 2005). One 
noticeable perturbation in a rotation curve occurs when the 
velocity on one side decreases considerably while the other 
side remains fairly constant. Commonly these features, so 
called bifurcations, happen when there is a close companion 
galaxy. An example of this behaviour can be seen in Fig. 
[T] obtained by Fuentes Carrera et al. (2004), where the RC 
of a galaxy in a pair of galaxies of similar sizes and masses 
is displayed. 
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The natural question to ask here is how these velocity 
disturbances relate to the passage of the companion galaxy 
and if they can be used as a timer of the stage of evo lution 
of an interaction. Recently, Kronberger et al. (2006) (here- 
after Kron06) have investigated disturbed RCs in siniulated 
interacting galaxies. Using simulations, those authors re- 
produced the observations of intermediate redshift galax- 
ies {z^O.5) by putting a virtual slit along the major axis 
of each investigated system (see also Barton et al 1999|) . 
They also found that the features associated with the per- 
turbations in the RCs strongly depend on the viewing an- 
gles, particularly, for the derived asymmetry parameters 
([Dale et al. 200ip . which presented quite a large disper- 
sion, for minor mergers and fly byes. For major mergers, no 
similar detailed study of the temporal behaviour of these 
asymmetries has been made. This motivated us to revisit 
Kronberger's study focusing on the analysis of the time evo- 
lution of the asymmetries for major mergers. Furthermore, 
we chose to peforme our analysis in 3D, to uncover the real 
level of disturbance of the systems. For this purpose, we 
used N-body/hydrodynamic simulations with cooling, star 
formation and supernova feedback. 

2. The Numerical Experiments 

We ran pre-prepared initial conditions of major mergers 
involving Milky- Way type galaxies. We used the proce- 
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Fig. 1. Rotation curve of the spiral galaxy NGC 
5953 in the interacting galaxy pair KPG 468 
(| Hernandez Toledo et al. 2003P within the inner 15", 
derived from a full 2D velocity map. A velocity bifurcation 
between the approaching side (full circles) and receding 
one (empty circles) can be seen at 9". The companion 
NGC 5954 is a galaxy of similar size and morphological 
type 



dure developed by Springel et al. (2005) to construct near- 
equilibrium galaxies, each one consisting of a dark matter 
halo, a disk of gas and stars and a bulge, with a virial ve- 
locity of 160 km s~^. The gaseous and stellar disks have 
exponential surface densities with an scale- length, hd, of 
2.5 kpc h~^ and scale-height of 0.2 hd- The bulge com- 
ponent follows a Hernquist profile ( [Hernquist 1990p with 
scale-length of 0.2 hd- The dark matter halo is consistent 
with a Hernquist profile scaled to match the inner density 
distribution of dark matter haloes found in cosmological 
simulations (jNavarro et al. 1996P with a concentration pa- 
rameter c = 9. We use 20000, 10000, 50000 and 30000 par- 
ticles to represent the stellar disk and bulge, the gas and 
the dark matter components, respectively, of a given galaxy. 
Baryonic masses range from 1.9 x 10^h~^ Mq for particles 
on the disk to 1.3 x 10^h~^ Mq in the bulge component, 
while the dark matter particles have 3 x lO^h"-"^ -^0- We as- 
sume a 10% fraction of baryons initially in the form of gas. 
A higher resolution simulation with 600000 total particles 
was also run to discard resolution effects. 

The two identical galaxies were given an initial rel- 
ative velocity of 200 km s^^. The minimum separation, 
Tp, achieved after 1.3 Gyr h~^, during the first pas- 
sage is 30 kpc h^^. Galaxies are set on parabolic orbits 
and will merge due to dynamical friction. We investi- 
gated four possible configurations for the encounter. We 
ordered the experiments us ing the orbital parameters of 

as ii = and 12 = for A, 

«i = 45, 12 = -45 and uj2 = 90 

D. These ini- 



Toomre and Toomre (1972) 
ii = and ^2 = 180 for B 
for C, ii = -45, 12 = 45 and 0J2 = 270 for 
tial configurations were chosen to study the kinematical 
behaviour of the galaxies between the extreme cases of a 
parallel encounter versus a perpendicular one. For experi- 
ment A we also run two simulations changing the minimum 
separation at pericentre: Al and A2 with r^ = 15 kpc h^^ 
and Tp = 50 kpc h~^, respectively. 

The simulations were run by using the Tree-SPH 
Gadget-2 code (Spri ngel 2005[ ). The cold gas can form new 
stars following the Kennicutt law according to a stochas- 
tic algorithm (Scannapieco et al. 2005). However, since the 
initial gas component is a small fraction of the total mass. 



the impact of new stars is not important for this analy- 
sis. The systems experienced a number of close encounters 
until they finally merged. In this work we focus only on the 
perturbation arising during the first close passage. 



3. Simulated Rotation Curves 

In a disk structure supported by rotation and in equilib- 
rium within its potential well, the only velocity component 
should be a tangential velocity {v^) on the plane perpendic- 
ular to the angular momentum of the system. In practice, 
this is not forcibly the case. Motions of stars and gas in the 
disk can be modified by the presence of structures in the 
disk or by external perturbations. In this case, particular 
features on the RCs, such as bifurcations, may appear. 

At a given time during the interaction, the simulated 
RCs are constructed by calculating v^ on the rotation plane 
(defined as the plane perpendicular to the total angular mo- 
mentum of the system) as a function of the distance to the 
mass centre of the galaxy (r). For the purposes of analysing 
possible asymmetries in the RCs, we considered v^ within 
four cones of solid angle of 60 degrees along two orthonor- 
mal directions on the disk plane: a — /?, 5 — ^. We note 
that results are insensitive to variations of these angles. In 
fact, we have varied these angles between 10 and 60 degrees, 
finding no significant changes in the trends. By analysing 
the velocities and angular momentum of the mass within 
these angles, we can study the effects of interactions in dif- 
ferent regions of the disks and look for loss of symmetry in 
the velocity patterns. 

We calculated the mean tangential velocity, <v^>, es- 
timated in equally spaced radial bins of size 1.5 kpc /i~^. 
We also estimated the predicted velocity that particles on 
circular orbits within a potential well dominated by dark 
matter should statistically have as Vcir{r) = [2Misil<llf.b 
where Mtot(< t) is the total mass (baryons plus dark mat- 
ter) within r. For systems in rotational equilibrium within 
such a potential well, we should find that < z;^ >~ V^ir. 
Important departures from this relation will be assumed to 
be indicative of a loss or a gain of angular momentum by 
the action of an external force. Before the galaxies start 
interacting, the disks are in rotational equilibrium within 
their potential wells, regardless of the angular sector chosen 
to measure <v^>. 

In order to assess the effects of the interaction on the 
rotation curves, we analysed the departures of < f ^ > from 
V^ir within each adopted angle as a function of time, for the 
four experiments performed. Fig. [2] shows an example of a 
disturbed RC in the a — f3 direction for experiment A. A 
very asymmetric response of the disk can be seen. For this 
RC, the velocity on the left side of the galaxy (/? direction) 
decreases considerably, while the velocity on the right-side 
of the galaxy (a direction) remains fairly constant, produc- 
ing a bifurcation similar to those reported by observations 

(Fig. m- 

In order to study the presence of these bifurcations in 
the RCs along the interaction event, we defined the resid- 
ual velocities between the RC and the circular velocity as: 
R =< W0 > — Vcir. By estimating i? as a function of radius, 
we can quantify the difference - and its amplitude - between 
the measured tangential velocity on one side of the galaxy 
and that on the opposite side. 
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Fig. 2. Circular (continuous line) and tangential velocities 
(triangles and squares) in experiment A showing clear bifur- 
cation. Triangles and squares denote opposite orientations 
along a given direction used to measure the curves. The side 
which shows a declining tangential velocity is the one clos- 
est to the companion galaxy. Error bars correspond to the 
standard deviations of the mean values computed within 
each bin. 
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Fig. 3. Velocity residuals as a function of time estimated 
for baryons within a (squares) and /? (triangles) angles for 
experiment A. 



Since we are interested in the evolution of these differ- 
ences with time, we defined a mean residual velocity < R > 
at each given time, averaging out the corresponding residu- 
als (the central regions of the galaxies are excluded to avoid 
an artificial weakening of the signal originated by the disk 
structure). Fig. [3] shows the mean residuals calculated for 
experiment A as a function of time in the a — f3 direction. 
From this figure, we can observe that, in general, significant 
residuals are present, starting at t « 1.25 Gyr h~^, which 
corresponds to the time of peripassage. We note that the 
largest bifurcations appear when the RC is measured within 
the angular sector closer to the companion galaxy. This can 
be seen in Fig. [3] where the /3 direction, which is the one 
closest to the companion, shows the largest mean residu- 
als. A detailed analysis of the angular momentum content 
of the simulated galaxies shows that these residuals can be 
directly linked to the variation of angular momentum on 
orthogonal directions. 



4. Can bifurcations be used as a timer for recent 
interactions? 

In order to quantify if bifurcations in the RCs can be used 
as a timer of the stage in the interaction, we analyse the 
evolution with time of the global RC asymmetry, defined 
by Dale et al. (2001). This asymmetry parameter (here- 
after AP) is calculated as the area between the approach- 
ing and receding halves of the RC, normalized by the av- 
erage area under the curve. Fig. 3] shows that when the 
simulated galaxies have galaxy-galaxy distances of approx- 
imately r = 30 kpc h~^, significant signal of asymmetry 
appears. Similar levels of AP are detected in our four ex- 
periments, although the details of the features vary from 
one to the other. 

Hence, according to the trends shown in Fig. [3] and 
Fig. m two galaxies need to be very close to start experi- 
encing a significant, internal re-distribution of angular mo- 
mentum and mass, but when this has occurred, the RCs 
remain perturbed for an extended period until they reach 
equilibrium again, after almost 1 Gyr h~^ of the pericen- 
tre, when the centres of masses are separated by more than 
200 kpch^^, in agreement with previous findings for minor 
mergers and fly- byes (e.g. Kron06). 

A detailed analysis of the evolution of the AP displayed 
in Fig. m shows that before the peripassage, for the four 
experiments remains approximately at a value of 5%. From 
the peripassage, a noticeable increase is detected with val- 
ues up w 20% and greater than 10%, and with a mean AP 
that, at least, doubles the level registered up to peripassage. 
After 0.5 Gyr h^^ from peripassage, it starts decreasing. A 
similar analysis for Al and A2 shows that AP anticorrelates 
with Tp. We find no AP greater than 10% for Vp larger than 
~ 50 kpc h^^ in this merger configuration. 

We calculated the AP for the RCs of 3 galaxies in iso- 
lated pairs of galaxies of similar masses, studied by Fucntes- 
Carrera (j Hernandez Toledo et al. 2003^ and 2 galaxies in 
Hickson Compact Groups observed by (Mendes de Oliveira 
et al. 2003) which show bifurcations. All RCs were de- 
rived from full 2D velocity field maps by averaging large 
velocity cones around the major axes of the galaxies, and 
therefore they are good representations of their kinemat- 
ics as a whole. They have been corrected by inclination 
so that the comparison with simulations is as good as it 
can be carried out. In the case of the pairs studied by 
Fuentes Carrera the asymmetry values obtained are: 19% 
for NGC5426, 13% NGC5427 and 9% for NGC5953. A sim- 
ilar behaviour is obtained for the galaxies in groups: 13% 
for HCG7c and 36% for HCG96a. We have also calculated 
AP for non-bifurcated-RC galaxies (HCG lOd, HCG88a, 
HCG 88c, HCG 88d, HCG 89a, HCG 91cl, HCG 96d) in 
compact groups. The obtained mean value was 7%, which 
is very similar, qualitatively, to the one found in the simula- 
tions, during the unperturbed phase (5%). It is interesting 
to note that the galaxies with bifurcated RC in compact 
groups are spiral galaxies with only one companion of sim- 
ilar mass (HCG 07a and HCG 96c, respectively), and sev- 
eral other considerably less massive companions (1/3 of the 
mass of the first ranked galaxy, at most). 

Hence, from the analysis of these high quality observed 
RCs we get similar results to those found in simulations for 
the asymmetry value of the bifurcated RCs. 
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Fig. 4. Asymmetry parameter as a function of time for our four experiments (A, B, C, D) along the a — (3 directions 
(see Section 2 for experiment description). 



5. Conclusions 

We studied the RCs of two simulated Milky Way type 
galaxies during their first encounter. We performed four 
simulations with different orientations of the angular mo- 
mentum direction of the systems. The simulated RCs dis- 
played bifurcations consistently with observational findings. 
The evolution with time of the features were investigated. 
Our main conclusions are: 



1 



The asymmetry parameter (defined bv lDale et al. 200T|) 
shows a clear increase at the peripassage, reaching val- 
ues larger than a factor of two the levels detected for 
RCs in equilibrium. After w 0.5 Gyr h~^, the asymme- 
try parameter begins to decrease again to lower values. 
We found that observed RCs corrected by inclination 
yield a similar behaviour. 
2. The most significant asymmetries are detected at peri- 
passage when the systems are separated ~ 30 kpc h~^. 
In our experiments we found that the perturbations in 
the RCs can be detected even after almost 1 Gyr h~^ 
of the peripassage in agreement with previous works for 
minor mergers and fly byes (Kron06) 

Therefore the asymmetry parameter could be used to 
roughly estimate if the pericentre has occurred in the last 
f» 0.5 Gyr h^^. The question if these features can always 
be seen in RCs of interacting galaxies might depend on 
the viewing angle (Kron06). However, the good agreement 
of the observational and numerical results obtained in this 
paper encourages us to carry out a more comprehensive 
comparison with observations in a future paper. 
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